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Abstract 

Yukawa couplings between Higgs bosons and quarks in the Standard Model (SM) 
and supersymmetric theories can be measured in the processes e + e~ — > QQ + Higgs. 
The cross sections and Higgs energy distributions of these processes in the SM and 
minimal supersymmetric model have been determined including the complete set of 
next-to-leading order QCD corrections for all channels. 
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1 Introduction 



In the Standard Model (SM) and its supersymmetric extensions, the masses of electroweak 
gauge bosons, leptons, and quarks are generated by interactions with Higgs fields The 
Yukawa couplings between Higgs particles and fermions therefore grow with the masses Mf 
of the fermions. The couplings obey a universal scaling law g/fH = Mf/v in the SM, with 
v pa 246 GeV being the ground-state value of the Higgs field. In supersymmetric models, 
which involve at least two Higgs doublets, the size of the Yukawa couplings is also set 
by the fermion masses, yet the relationship is more complex owing to the mixing among 
the Higgs fields. The Yukawa couplings @ of the two CP-even light/heavy Higgs bosons 
h/H and of the CP-odd Higgs boson A in the minimal supersymmetric extension of the 
Standard Model (MSSM) || Q, expressed in units of the SM couplings, are determined by 
the parameters tg/3 = v 2 /vi, the ratio of the vacuum expectation values of the two Higgs 
fields generating the masses of up- and down-type particles, and a, the mixing angle in 
the CP-even sector. In the decoupling limit, in which the light Higgs mass reaches the 
maximum value for a given parameter tg/3, the h Yukawa couplings approach the SM 
values. In general, they are suppressed for up-type fermions and enhanced for down-type 
fermions; the enhancement increases with tg/3 and can therefore be very strong. 
Higgs radiation off top or off bottom quarks in e + e~ collisions, 

e+e" _> QQcf) [Q = t,b; <p = H SM , h, H, A], (1) 

lends itself as a suitable process for measuring the Yukawa couplings of all Higgs bosons 
||, particularly for the SM Higgs bosons, the light Higgs boson h and for moderately heavy 
Higgs bosons H and A. Within the SM it has been demonstrated that the top Yukawa 
coupling can be measured with an accuracy of 2-3% at linear colliders with yfs = 800 
GeV ||. We present the cross sections for these processes including the next-to- leading 
order (NLO) QCD corrections [[?], ||. In this analysis we have calculated the complete set 
of 0(a s ) QCD corrections to all subchannels in ([I]) systematically. The large number of 
interfering subchannels in supersymmetric theories renders this program more complex 
than the corresponding calculation in the SM, in particular since the relative weight of 
the subchannels varies over the supersymmetric parameter space and over the phase space 
for different mass ratios. Moreover, the Higgs-energy distributions haven been obtained 
at NLO |J. They turn out to be relevant in the separation of the resonant parts and 
those, which depend on the Yukawa couplings. Introducing appropriate cuts in the Higgs 
energy the sensitivity to the Yukawa couplings can be increased experimentally. 

2 QCD Corrections 

The QCD corrections can be categorized into five classes. Virtual corrections of the 
internal quark lines, of the 7/Z-quark vertices, of the Higgs-quark vertices, and box 
diagrams interfere with the Born amplitude. Gluon radiation off internal and external 
quark lines adds incoherently to the cross sections. The value of the electromagnetic 



1 



coupling is taken to be a = 1/128 and the Weinberg angle as sin 2 6w = 0.23. The mass 



of the Z boson is set to M z 
quarks are set to M t = 174 GeV 



91.187 GeV, and the pole masses of the top and bottom 
T0| and[] Mb = 4.62 GeV Jll|, respectively. The masses 



of the MSSM Higgs bosons and their couplings are related to tg/3 and the pseudoscalar 
Higgs boson mass Ma- In the relation we use, higher-order corrections up to two loops 
in the effective-potential approach are included [12|. The SUSY parameters are chosen as 
[i = A t = Ab = and Mq = 1 TeV; this simple choice is sufficient to illustrate the main 
results. 

The production of bb<fi final states can be mediated by resonance channels e + e _ — > 
Ism 



ZH SM} Zh, ZH and e + e~ — > Ah, AH. We describe the resonance structures as Breit- 
Wigner forms by substituting M 2 — > M 2 — iMT in all boson propagators. The decay 
widths of the Higgs bosons are calculated including higher-order corrections, as described 
in Refs. fl3| , 14}] , while the Z width is set to Tz = 2.49 GeV. For tt(f) production the 
widths can be neglected, since the Higgs masses are taken below the ti threshold. The 
renormalization scale of the QCD coupling a s , which is evaluated in NLO with five active 
flavors normalized to a s (M z ) = 0.119 |TI|, is chosen at /j, 2 r 



s, where s 



E% u is the 



0.119 [ 

center-of-mass (CM) energy squared. 

The QCD radiative corrections have been calculated in the standard way. The Feyn- 
man diagrams have been evaluated within dimensional regularization. Ultraviolet diver- 
gences are consistently regularized in D = 4 — 2e dimensions, with 75 treated naively since 
no anomalies are involved. The renormalization of the QQ<p vertices is connected to the 
renormalization of the quark masses, which, in the case of the top quark, is defined on 
shell (see e.g. Refs. [15, 16| for details). In the case of bottom quarks large logarithms are 
mapped into the running mass Tn^Qmggs) for the Yukawa couplings of the b quark, with 
Qffiggs denoting the squared momentum flow through the corresponding Higgs-boson line. 
The infrared divergences encountered in the virtual corrections and in the cross section for 
real gluon emission, are treated in two different ways. Both calculations follow subtrac- 
tion procedures, one using dimensional regularization and one introducing an infinitesimal 
gluon mass |17|. The results obtained by the two different procedures are in mutual nu- 
merical agreement after adding the contributions from virtual gluon exchange and real 
gluon emission. A second, completely independent calculation of the QCD corrections 
to the total cross section was based on the evaluation of all relevant cut diagrams of the 
photon and Z-boson self-energies in two-loop order, generalizing the method applied to 
tt(g) intermediate states in Ref. [18]. The results of the two approaches are in numerical 
agreement. Moreover, the parts that have been calculated in Refs. [[L^, |2(|, i.e. photon 
exchange and resonant contributions, are in full agreement with our corresponding partial 
results. 



1 This value for the perturbative pole mass of the bottom quark corresponds in NLO to an MS mass 
m h {m b ) = 4.28 GeV. 
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3 Results 



a.) Asymptotic behaviour. The QCD corrections to the top final states ticj) can be 
interpreted easily in two kinematical areas. Whenever the invariant mass of the tt pair is 
close to threshold, the gluonic Sommerfeld rescattering-correction is positive and becomes 
large. In the threshold region the K factor approaches the asymptotic form [7| 



K, 



tt4> 

thr 



1 + 



32a s 
Wt 



(2) 



with the maximal quark velocity fit = y (v 7 ^ — ^>) 2 — 4M 2 /2M 4 in the (tt) rest frame. 

For high energies, on the other hand, the QCD corrections are of order a s /ir. In the 
energy region s ^> 4M 2 3> M\ but log s/Mf ^> 0(1), which is relevant for the present 
analysis, the QCD corrections can qualitatively be traced back to vertex corrections and 
infrared gluon radiation. Since scalar Higgs bosons are radiated off top quarks preferen- 
tially with small energy [x — E^,jE t — > 0], as is evident from the leading (universal) part 
of the fragmentation function 



f(t^tH;x) 



9ttH 
16/T 2 



1 — X s 

+ X l °STj2 



(3) 



the QCD correction of the scalar Yukawa vertex, regularized by soft gluon radiation, 
approaches the value []19 | 
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The scalar Yukawa vertex is therefore reduced by four units in a s /7r which are compen- 
sated only partly by one unit due to the increase of the tt production probability, leading 
in total % 011 to 

Kit - 1-3— + ... ix.ro = H.u.h.H . (5) 



1-3— + 

7T 



The ellipsis accounts for hard Higgs and gluon radiation (of order +a s /Tr). Thus, the QCD 
corrections are expected negative for scalar Higgs particles in the high energy continuum. 

By contrast, the corresponding fragmentation function for the pseudoscalar Higgs 
boson |?TJ 

s 



f{t^tA-x) 



i*k x log 

16tt 2 g M 2 



(6) 



is hard so that the average of the vertex and IR gluon corrections over the Higgs spectrum 
amounts to 
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(7) 



Adding to this correction the increase of the tt production probability of one unit, the K 
factor is very close to unity 
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ttA 
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After hard gluon bremsstrahlung is taken into account (symbolized by the ellipsis), the 
overall QCD corrections for the pseudoscalar Higgs boson are therefore expected slightly 
positive. [For ultra-high energies, i.e. log s/Mf 3> 1, hard gluon bremsstrahlung becomes 
important. Similarly to the leading terms in the fragmentation functions eqs. (j3|) and 
the QCD corrections for scalar and pseudoscalar Higgs bosons approach each other as a 
result of chiral symmetry restoration in asymptotia; this has been verified in a numerical 
calculation.] 

Similar estimates can be applied to bottom final states which in general are dominated 
by resonance decays. After absorbing the large logarithms loglQ^/M^) into the Yukawa 
couplings, the non-leading effects are positive: 

iY r 6 ^«l + {¥,l}^ for {Higgs, Z}^ bb . (9) 

Also close to the thresholds and in the high-energy limit the QCD corrections remain 
positive after mapping the large (negative) corrections into the running Yukawa couplings. 
Since different channels are activated at the same time, only a qualitative estimate can 
be given in the continuum regime, 

Kti = l + c^ withe = 0(1), (10) 

7T 

while details must be left to the numerical analysis. 

b.) Numerical results. The total cross sections for Hti production in the SM are 
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Figure 1: The total cross section for the process e + e _ — > ttH+X , including QCD radiative 
corrections (full curve) and at LO (dashed curve) as a function of the scaled Higgs energy 
x H /§/. 



presented in Fig. [I] as a function of the Higgs mass for different collider energies M. The 
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NLO (LO) cross sections are given by the full (dotted) curves. The QCD corrections are 
large for y/s = 500 GeV due to the Coulomb singularity at threshold and moderate for 
high energies. At high energies the results agree with the approximation of eq. (|5]). The 
cross sections amount to more than about 0.1 fb, which leads to a significant number of 
events at the TESLA collider, being designed to reach integrated luminosities of about 
/ C ~ 1 ab _1 in three years of operation. The corresponding Higgs-energy distribution is 
presented in Fig. |2] for M H = 120 GeV as a function of the scaled variable xh = 2E H / y/s 
||. The shape of the distribution is shifted towards larger values of xh due to the Coulomb 
singularity at threshold. The results for small values of xh reproduce the approximation 
ofeq. @. 
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Figure 2: The Higgs-energy distribution for the process e + e~ — > tiH + X , including QCD 
radiative corrections (full curve) and at LO (dashed curve) as a function of the scaled 
Higgs energy x H 

The results for the MSSM are exemplified in Fig. ||] for tg/3 = 30 and 30 and for the 
collider energy Ecu = 500 GeV [§] . If required by the size of the cross section, which 
should not fall below ~ 10 -2 fb in order to be accessible experimentally, we switched to 
Eqm = 1 TeV. The Born terms are shown by the dotted curves, while the final results for 
the cross sections, including QCD corrections, are given by the full curves. 

For E C m = 500 GeV the QCD corrections to Higgs-boson production in association 
with ti pairs increase the scalar Higgs-production cross sections significantly, as can be 
inferred from Fig. || Close to threshold the numerical results clearly exhibit the strong 
increase of the cross sections due to the Coulomb singularity (0). Moreover, for tg/5 = 30 
the cross sections are strongly suppressed except for the regions where the light (heavy) 

2 Although the value tg/3 = 3 is already excluded for vanishing mixing in the stop sector, this choice 
exemplifies the complete result for smaller values of tg/3 and can easily be extended to the non- vanishing 
stop mixing case. 
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Figure 3: Production cross sections for the MSSM Higgs bosons h, H and A in association 
with heavy t, b quark pairs fi^j; Born approximation: dashed, QCD corrections included: 
full curves. The rapid drops in the bbA cross section at tg/3 = 3 are due to the kinematical 
opening of the resonant H — > WW, hh decays, which reduce the branching ratio of the 
resonant H — » bb decay. 
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scalar Higgs mass is close to its upper (lower) bound. For tg/3 = 3 the cross section 
amounts to about 1 fb, which leads to a significant number of events at the TESLA 
collider. 

For CM energies of 1 TeV the QCD corrections to scalar Higgs production in associa- 
tion with ti pairs are of moderate size. They decrease the cross sections by about 3-5%. 
This is in accordance with the asymptotic form of the K factor . For pseudoscalar Higgs 
production, the size of the QCD corrections is slightly positive at 1 TeV, in agreement 
with the qualitative discussion above. 

The cross sections for Higgs-boson production associated with bb pairs are significantly 
larger due to the resonance contributions from on-shell Z and Higgs-boson decays into 
bb pairs. The QCD corrections increase these cross sections by about 5-25%. The drop 
in the bbA production cross section for tg/3 = 3 at Ma ~ 140 GeV and 175 GeV can be 
attributed to the crossing of the thresholds for resonant H — > WW and H — > hh decays, 
respectively, in HA final states. 

Without cuts, the intermediate resonance decays Z, h,H,A — > bb dominate all bb<p 
production processes, whenever they are kinematically allowed, and it will be difficult 
to extract the bottom Yukawa couplings experimentally in regions where resonant Higgs 
decays to bb pairs are dominant. This is the case at large values of tg/3 for all neutral 
Higgs particles and at small values of tg/3 for Higgs masses below the WW threshold. In 
these cases the branching ratios, which determine the size of the bb<p cross sections, will 
be nearly independent of the bottom Yukawa couplings. Resonance decays R — *> bb in the 
bbh/H/A final states can however be eliminated by cutting out the resonance energy of the 
final-state Higgs boson, £^ >res — (s + M| — M|,)/2- v /i. After subtracting these resonance 
parts, the non-resonant contributions are suppressed by about one to three orders of 
magnitude. The resonances have been removed from the cross sections in the examples of 
Fig. H by subtracting the two-particle cross sections in the Breit-Wigner bands Mr ± A 
of the energy with the resolution A = 5 GeV. This theoretical definition is used for 
the sake of simplicity; wider cuts may be required in experimental analyses. Peaks and 
dips in the cross sections are the result of overlapping Breit-Wigner bands. For the scalar 
Higgs particles they arise from overlapping Z and A boson bands; in pseudoscalar Higgs 
production the light and heavy scalar resonance bands overlap for 100 GeV < Ma < 120 
GeV for tg/3 = 30. The dips occur whenever the two resonance bands touch each other, 
while the peaks between the dips occur when the resonance masses coincide exactly. As 
shown in Fig. ^, the QCD-corrected cross sections are still close to 1 fb or slightly below, 
except for heavy masses at small values of tg/3. Thus, ensembles of order 10 3 events can 
be collected at a high-luminosity collider |5|. 

Fig. H presents the pseudoscalar Higgs-energy distribution in tiA production at y/s = 1 
TeV for Ma = 120 GeV and tg/3 = 3 as a function of the scaled Higgs energy xa = 2Ea/ y/s 
0. As in the case of the SM Higgs boson the distribution is shifted towards larger Higgs 
energies due to the Coulomb singularity at threshold. Moreover, it is clearly visible from 
the comparison of Figs. |2] and |^ that the pseudoscalar energy spectrum is harder than the 
scalar one, which is related to the absence of the 1/x enhancement in eq. (||). 

The corresponding Higgs-energy distributions of bbh production for Mh = 110 GeV is 
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Figure 4: Continuum production of MSSM Higgs bosons h, H, A in association with a bb 
pair after subtraction of resonance decays to bb pairs in the Breit- Wigner bands Mr ± 5 
GeV®. 
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Figure 5: The Higgs-energy distribution for the process e + e~ — > tiA + X , including QCD 
radiative corrections (full curve) and at LO (dashed curve) as a function of the scaled 
Higgs energy xa for Ma = 120 GeV at yfs = 1 TeV ||/. 
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exemplified in Fig. || and for bbA production for Ma = 120 GeV in Fig. |^ as functions of 
the corresponding scaled Higgs energies for tg/3 = 30 at yfs = 500 GeV . The resonance 
peaks of A, Z — > bb are dominating the distribution of bbh production and of h, H — > bb the 
distribution of bbA production. The continuum contributions are suppressed by several 
orders of magnitude. The QCD corrections are moderate for all Higgs energies. 
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Figure 6: The Higgs-energy distribution for the process e + e — > bbh + X , including QCD 
radiative corrections (full curve) and at LO (dashed curve) as a function of the scaled 
Higgs energy x h for M h — 110 GeV at sfs = 500 GeV $J. 



4 Conclusions 

Top- and bottom- Yukawa couplings in the SM and supersymmetric theories can be mea- 
sured at future linear e + e~ colliders in Higgs radiation off top and bottom quarks. We 
have presented the total cross sections and Higgs-energy distributions including the full 
NLO QCD corrections. The corrections turn out to be large for Higgs radiation off top 
quarks at cm. energies y/s = 500 GeV, while they are moderate at larger energies. The 
QCD corrections to Higgs radiation off bottom quarks are moderate after absorbing large 
logarithms in the running bottom Yukawa couplings. 

Measurements of the top Yukawa coupling in SM Higgs radiation off top quarks can be 
performed quite accurately at the future TESLA collider thanks to the high luminosities. 
Measurements of Yukawa couplings in supersymmetric Higgs radiation off heavy quarks 
at e + e~ linear colliders are difficult. This is a result of the large number of subchannels 
contributing to the QQh/H/A final states in supersymmetric theories in general, and the 
contamination by two-Higgs final states in particular. Nevertheless, the continuum cross 
sections appear large enough to allow for a solution of this experimental problem as shown 
in the present analysis. Even though experimental simulations are beyond the scope of 
this note, it may be concluded from earlier SM analyses that the method will work at 
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Figure 7: The Higgs- energy distribution for the process e + e — ► bbA + X , including QCD 
radiative corrections (full curve) and at LO (dashed curve) as a function of the scaled 
Higgs energy x A for M A = 120 GeV at = 500 GeV gj. 



least in part of the MSSM parameter space, thus providing us with the absolute size of 
the quark-Higgs Yukawa couplings in the minimal supersymmetric theory. 
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